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We have applied our recent approach (Kargarian, et.al Phys. Rev. A 76, 60304 (R) (2007)) 
to study the quantum information properties of the anisotropic s=l/2 Heisenberg chain. We have 
investigated the underlying quantum information properties like the evolution of concurrence, en- 
tanglement entropy, nonanalytic behaviours and the scaling close to the quantum critical point of 
the model. Both the concurrence and the entanglement entropy develop two saturated values after 
enough iterations of the renormalization of coupling constants. This values are associated with the 
two different phases, i.e Neel and spin liquid phases. The nonanalytic behaviour comes from the 
divergence of the first derivative of both measures of entanglement as the size of system becomes 
large. The renormalization scheme demonstrates how the minimum value of the first derivative and 
its position scales with an exponent of the system size. It is shown that this exponent is directly 
related to the critical properties of the model, i.e. the exponent governing the divergence of the 
correlation length close to the quantum critical point. We also use a renormalization method based 
on the quantum group concept in order to get more insight about the critical properties of the model 
and the renormalization of entanglement. 
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I. INTRODUCTION 

The quantum information theory (QIT) and condensed 
matter physics share in using the entanglement--. The 
main motivation behind such interest is two folds: On 
the one hand in QIT the entanglement is recognized as 
an essential resource for implementing some quantum in- 
formation tasks such as quantum computations, quantum 
cryptography and densecoding^ . It is also believed that 
the protocol based on the quantum entangled state has 
much speed than the classical ones. On the other hand 
entanglement is a unique measure of the quantum corre- 
lation of a pure state in condensed matter physics. Thus, 
for condensed matter systems the entanglement can de 
considered as some kind of correlation which appears in 
the pure states. 

The role of entanglement in quantum phase transi- 
tion (QPT)^ is of considerable interest-. Quantum phase 
transitions occur at absolute zero temperature. They are 
induced by the change of an external parameter or cou- 
pling constant, and are driven by quantum fluctuations. 
Quantum as well as classical phase transitions are charac- 
terized by detecting nonanalytic behaviors in some phys- 
ical properties of the system. It is often accompanied by 
a divergence in some correlation functions, but quantum 
systems possess additional correlations which do not ex- 
ist in a classical counterparts, the entanglement. Entan- 
glement is a direct measure of quantum correlations and 
shows nonanalytic behavior such as discontinuity in the 
vicinity of the quantum critical point^'®. In the past few 
years the subject of several activities were to investigate 
the behavior of entanglement in the vicinity of quantum 
critical point for different spin model3^'^'^ i^'^°i^^ as well 
as itinerant system o^^'^'^i^^ . 

In our recent worki^ we have introduced the notion of 



renormalization of entanglement and showed that this 
notion truly describe the nonanalytic behavior of the 
derivative of entanglement close to critical point of the 
Ising model in Transverse Field (ITF). Moreover, we have 
investigated how the critical point is reached by increas- 
ing the size of system via the renormalization group (RG) 
approach. The finite size scaling demonstrates that the 
RG of entanglement truly captures the critical behavior 
of the model. The renormalization of quantum states has 
also been introduced in terms of matrix product statesi^. 

In this work we have applied the same approach to the 
anisotropic spin 1/2 Heisenberg (XXZ) model in order to 
investigate the critical behavior of the model using the 
evolution of entanglement through the quantum renor- 
malization group^ii^ close to quantum critical point. 
In this model the concurrence between the two nearest- 
neighbor sites develops a maximum at the isotropic point 
(A = 1) without any singularity in its first derivatives^, 
i.e it vanishes at the critical point A = 1. However, in our 
approach a singularity, i.e the divergence of first deriva- 
tive, is seen at the critical point. It is expected that our 
approach is efficient and powerful to study the quantum 
phase transition in various systems since the thermody- 
namic limit of the model is simply captured through an 
analytic RG equations. Quantum renormalization group 
with fixed boundary conditions give rise only to a non- 
trivial fixed point A = 1 and trivial fixed points A = 0, cxd 
for the XXZ model while for every point in the < A < 1 
region the model is critical, i.e gapless. A modification 
of the boundary conditions by means of quantum group 
concept^° restores the gapless properties of < A < 1, 
although the critical behavior of the model and its scaling 
properties remain unchanged from the gapped phas. 

This article is organized as follows. In section-llll the 
main idea of quantum renormalization group is briefiy re- 
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viewed and in section- HTTl this idea is applied to the XXZ 
model. Sections- IIVI discusses the main idea of this paper 
about the renormalization of entanglement entropy and 
concurrence is introduced and the section-TVl is devoted 
to the quantum group analysis in order to get more in- 
sights about the critical features of the model. Finally in 
section- TVll we summerize our results. 



II. QUANTUM RENORMALIZATION GROUP 



The effective (renormalized) Hamiltonian is then 

H^f^ ^H^ff + Hf^ . (4) 

We will implement this approach in the next sections to 
obtain the quantum properties of the XXZ spin chain. 

III. RENORMALIZATION OF THE XXZ 
MODEL 



The main idea of the RG method is the mode elimina- 
tion or thinning of the degrees of freedom followed by an 
iteration which reduces the number of variables step by 
step until reaching a fixed point. In Kadanoff 's approach, 
the lattice is divided into blocks. Each block is treated 
independently to build the projection operator onto the 
lower energy subspace. The projection of the inter-block 
interaction is mapped to an effective Hamiltonian [H'^^^) 
which acts on the renormalized subspace^ii^.^ . The pro- 
cedure starts by decomposing the Hamiltonian into two 
parts, 



H 



BB 



(1) 



where the block Hamiltonian is a sum of commut- 
ing terms, each acting on different blocks. The inter- 
block interaction is represented by H^^ . The perturba- 
tive implementation of this method has been discussed 
in Refs. (21it24<] . We will shortly present this approach 
in the 1st order correction. The zeroth order effective 
Hamiltonian is given by 



H^" = PoH^Po, 



(2) 



where Pq is a projection operator. Since is a sum of 
disconnected block Hamiltonians 



N' 



1=1 



one can search for a solution of Pq in a factorised form 



N' 



^0 = n ^0 , 



1=1 



where N' is the number of blocks. In the standard quan- 
tum renormalization group, Pq is given by 



4 = EiV'«)(^'i' 

i=l 

where \'ipi) (i — 1, • • ■ ,k) are the k lowest energy eigen- 
states of hf . The interaction between blocks define the 
1st order correction by the following equation 



(3) 



The Hamiltonian of XXZ model on a periodic chain of 
N sites is 



H(J,A) 



7 ^ 



i+1 



A<a^+i) (5) 



where J, A > 0, J is the exchange coupling, A is the 
anisotropy parameter and cr" {a = x, yz) are Pauli ma- 
trices. For A = 1, the Hamitonian is SU (2) symmetry in- 
variant, but for A 7^ 1 the SU(2) symmetry breaks down 
to the U{1) rotational symmetry around the z-axis. The 
model is exactly solvable by means of Bethe Ansats as 
far as the rotational symmetry exists. It is known that 
for < A < 1 the model is gapless with qasi-long range 
ordered where the correlations decay algebraic with no 
magnetic long range order. For A > 1, the symmetrty 
is reduced to Z2 and a gap opens which is in the uni- 
versality class of one dimensional antiferromagnetic Ising 
chain. Indeed the third term in the Hamiltonian causes 
ordering in the system and as A tends to infinity the Neel 
state is the dominant phase of the system. The first two 
terms in the Hamiltonian extend the quantum fluctua- 
tions in the system and result in the corruption of the 
Neel ordering. It is shown that the competition between 
the quantum fluctuations and ordering yield a maximum 
value of concurrence between the two nearest neighbor 
sites which show a scaling behavior at the critical point 

^ _ ^25 26 

To implement the idea of QRG to calculate the en- 
tanglement and concurrence we use a three site block 
procedure where the block Hamiltonian is 



hf 



2.1 
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) + 



(6) 



In this case the inter-block (H^^) and intra- block (H^) 
Hamiltonians are 



N/3 



,(7) 



N/3 



(8) 
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where 



FIG. 1: Three sites in a block labeled sequentially. 



A remark is in order here, choosing the three site block is 
essential here to get a self similar Hamiltonian after each 
RG step. An odd site XXZ Hamiltonian has two de- 
generate ground states which are used to construct the 
projection operator of quantum RG (QRG). These de- 
generate ground states are 
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= (l TTi> + 9l TiT> + UTT)), 



(9) 



:(l Tii> + 9UTi> + UiT)), (10) 



where q is 



= -^[A + \/A2 + 8]. 



The corresponding energy is 



eo--j[A+VA2-f 8], 



(11) 



(12) 



and I T)) I i) are the eigenstates of . 

The projection operator (Pq) for the I-th block is de- 
fined 



Po =\t)ii{M + \i^)u{^ol 



(13) 



where | ff)/ and | JJ-)/ are renamed states of each block 
to represent the effective site degrees of freedom. The 
renormalization of Pauli matrices are given by 

Pq<iPo = (* - 1, 2, 3 ; a^x,y, z), (14) 

where 



2q 



fX,y _ 
?2 — 
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2 

2 + 
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(15) 



where the indices 1, 2 and 3 refer to sites labeled in a 
single block based on Fig.((T]). 

The effective Hamiltonian (using Eq.([l])) of the rcnor- 
malized chain is again an XXZ chain with the scaled cou- 
plings 



H-ff 



J' 



N/3 



,(16) 



A' = A- 



(17) 



The stable and unstable fixed points of the QRG equa- 
tions is obtained by solving A' = A = A*. The stable 
fixed points locate at zero and infinity while A* — 1 
stands for the unstable fixed point which specify the crit- 
ical point of the model. Starting with any initial value for 
A > 1, the coupling constant flows toward infinity show- 
ing that the model falls into the universality class of Ising 
model, while for A < 1 the stable A = fixed point is 
touched. As we have mentioned previously, for < A < 1 
the model represents a spin fluid phase. The transition 
between the two phases are truly captured by real space 
QRG. The main discrepancy of our results on the Ising 
model in transverse field^^ (ITF) and the present one on 
XXZ comes from the fact that XXZ model is critical for 
all values of < A < 1 but QRG equations do not to 
show the whole critical region except A = 0. In fact the 
QRG prescription only represents the masslessness prop- 
erty of the XY fixed point (A = 0). However, when the 
coarse graining procedure implemented by some appro- 
priate boundary conditions in order to get more corre- 
lation between blocks as is done in the quantum group 
method^*' , the critical line of the model is truly predicted 
(see section rv) . 



IV. RENORMALIZED ENTANGLEMENT AND 
ITS SCALING PROPERTY 

A. Entanglement and concurrence 

In this subsection we calculate the concurrence and en- 
tanglement of the XXZ chain using our proposal which 
implements the idea of renormalization group. As we 
have mentioned previously, a finite size block is treated 
exactly to calculate the physical quantities. The coupling 
constants of a finite block are renormalized via the QRG 
prescription to give the large size behaviour. The XXZ 
spin 1/2 Hamiltonian has two degenerate ground state 
in the three site block, however to define the density ma- 
trix we have to consider one of them. Thus, the density 
matrix is defined by 



= l'/'o)(0o|, 



(18) 



where {(po) has been introduced in Eq.®. The results 
will be the same if we consider |0q) to construct the den- 
sity matrix. 

There are basically two choices to define the concur- 
rence and entanglement for a three site block. {1) The 
symmetric case, in which the concurrence between site 1 
and 3 is obtained by summing over the degrees of free- 
dom of site 2. And the entanglement between site 2 and 
remaining sites of the block as measured by von-neuman 
entropy. (2) We sum over site 1 or 3 and get the con- 
currence between the middle site (2) and one site at the 
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FIG. 2: (color online) Representation of the evolution of con- 
currence in terms of RG iterations (steps). 



corner side of the block (see Fig.([T|)). Without loss of 
generality we only concentrate on case (1). The density 
matrix defined in Eq. (fT8|) is traced over site 2 degrees of 
freedom to get the reduced density matrix for site 1 and 
3 (91,3) which gives 



913 



/q2 0' 

110 

110 

V 0, 



(19) 



The corresponding eigenvalues oi R — 513^13 (where 
^13 — (^i ® o-|)^?*3(CTf ® cr|)) are in ascending order 



Ai — A2 — A3 — , A4 — 



Thus, the concurrence is obtained 



Ci3 = max{Ay^,0} 



(2 + g2)2 



(20) 



(21) 



According to Eg. (fTT|) . C13 is a function of A. The renor- 
malization of A defines the evolution of concurrence as 
the size of system becomes large. We have plotted in 
FIGlDthe value of C13 versus A for different QRG iter- 
ations. 

The plots of C13 versus A for different QRG steps 
cross each other at the scale invariant critical point, 
Ac = A* = 1. By increasing the size of system (higher 
QRG steps) the concurrence develops two different be- 
haviours which are separated at Ac. After enough QRG 
steps the value of C13 is saturated for < A < 1 which 
shows the existence of quantum correlations between two 
blocks at large distances in an infinite chain that is effec- 
tively described by a three site model with the renormal- 
ized coupling constants. In this region (0 < A < 1) the 



quantum fluctuations arising from the transverse inter- 
actions have dominant effect and destroy any long range 
order. In spite of lacking a long range order, the evolu- 
tion of concurrence via RG indicates that the spin fluid 
phase contains the quantum correlations, i.e the qubits in 
the presence of quantum fluctuations are quantum corre- 
lated. For A > 1, C13 saturates the zero value asymptot- 
ically to represent the lack of quantum correlation in the 
Ising limit. We have also examined the entanglement of 
formation^! between sites 1 and 3 which can be obtained 
by the following relation: 



El 



-yiog2(z/) - (1 - y) iog2(i - y), 



1 1 



y = 2 + 2^/^ 



^13- 



(22) 



The entanglement of formation shows similar behavior to 
FIGH] because of monotonous relation between concur- 
rence and entanglement of formation. In the symmet- 
ric case we can also study the entanglement between the 
middle site (2) with the remaining sites of the block. The 
amount of entanglement is given by von-Neuman entropy 
of reduced density matrix obtained after tracing out the 
remaining sites of block 



1 



Q2 = 



2 




(23) 



The von-Neuman entropy which is the entanglement of 
site 2 is 



2 , 2 

E = ^log. 



I -82 77—^ - log2 • (24) 



2 + "-^ 2 + g2 2 + g2 2 + g2 

The variation of entanglement (E) versus A has been 
plotted in FiglJl Different plots show the evolution of E 
under QRG iterations. In other words, the different step 
of QRG show how the entanglement evolves when the 
size of chain is increased. Similar to C13, ii^ behaves as an 
order parameter which gets a nonzero value for < A < 
1 and zero for A > 1 in the infinite size limit. Nonzero 
i? for < A < 1 verfies again that the state of model 
is entangled for < A < 1 where the ground state is 
characterized by a gapless excitation and algebraic decay 
of spin correlations (the spin fluid phase). While, E is 
zero for the Neel state (A > 1) which is not an entangled 
state. 

As previously discussed the entanglement between the 
middle site (2) and the remaining ones tends to maximal 
entanglement for the region of coupling < A < 1 after 
few RG steps FigO While the entanglement of forma- 
tion between sites 1 and 3 tends to a nonmaximal value, 
in fact this is a monogamy property of entanglement^^. 
On the basis of monogamy property , the maximal entan- 
glement between two parties restricts their entanglement 
with third party and vis versa. 
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FIG. 3: (color online) The evolution of Entanglement entropy 
as the size of system grows (i.e. different RG iterations). 

B. nonanalytic and scaling behavior 

The non-analytic behavior in some physical quantity is 
a feature of second-order quantum phase transition. It is 
also accompanied by a scaling behavior since the correla- 
tion length diverges and there is no characteristic length 
scale in the system at the critical point. There have been 
a number of theoretical studies on the entanglement and 
quantum phase transitions. These studies show that the 
entanglement as a direct measure of quantum correla- 
tions indicate critical behavior such as diverging of its 
derivative as the phase transition is crossed. Osterloh,et 
al^ have verified that the entanglement in the vicinity of 
critical point of ITF and XX model in transverse field 
show a scaling behavior. It has also been shown that the 
entanglement between a block of spins and the rest of sys- 
tem scales near the quantum critical point^. The scaling 
properties of entanglement with the size of block differs 
at and away from the critical point. At the critical point 
where the correlations decay algebraically entanglement 
scales logarithmically while it saturates away from the 
critical point^— . It can be interpreted in the framework 
of conformal field theory— associated with the quantum 
phase transition and the central charge of the theory. In- 
deed for the XXZ model, the entanglement entropy, i.e. 
the entanglement of a block (sublattice bipartition) with 
the rest of the system, demonstrates extermum behavior 
at the critical poini*^. 

In this work, we adopt a preferable and distinct way 
to study the block entanglement via the renormalization 
group approach. As we have stated in the RG approach 
for XXZ model, a large system, i.e. N = S""*"^, can be 
effectively described by three sites with the renormalized 
couplings of the n-th RG iteration. Thus, the entan- 
glement between the two renormalized sites represents 



the entanglement between two parts of the system each 
containing sites effectively. In this respect we can 
speak of block entanglement - block-block entanglement 
or the entanglement between a block and the rest of sys- 
tem. It is shown that the first derivative of the entangle- 
ment shows a diverging behavior as the critical point is 
reached. 

Having this in mind, the first derivative of concurrence 
is analyzed as a function of coupling A at different RG 
steps which manifest the size of system. The derivative of 
concurrence with respect to the coupling constant (^) 
shows a singular behavior at the critical point. It is given 

by 

dA ~ dA„ dA„_i "' dA ^ ' 

where stand for the concurrence or von-Neuman 

entropy and A„ the renormalized anisotropy coupling 
(A) at the n-th RG iteration. The singular behavior 
is the result of discontinuous change of C at A = A^,. 
We only concentrate on the nonanalytic behavior of en- 
tanglement between the middle site and remaining sites 
of each block. We have plotted ^ versus A in Figd] 
for different RG iterations which shows the singular be- 
havior as the size of system becomes large (higher RG 
steps). A more detailed analysis shows that the position 
of the minimum (A„i) of ^ tends towards the critical 
point like A^ = Ac -f- Af^^-*"^ which has been plotted in 
Figini Moreover, we have derived the scaling behavior 
of y = 1 A™ versus A^. This has been plotted in FiglH] 
which shows a linear behavior of ln(y) versus ln{N). The 
exponent for this behavior is | ^|a„ ~ N^''^'^ . This re- 
sults justify that the RG implementation of entanglement 
truly capture the critical behavior of the XXZ model at 
A = 1. 

Besides the entanglement entropy, two point quantum 
correlation as measured by concurrence develops a max- 
imum at the critical point. The concurrence of the two 
nearest neighbor sites of a large chain has numerically 
been calculated and the maximum at the critical point 
is due to the counterbalance of ordering and quantum 
fiuctuations^^. In the Ising limit, i.e A — > oo, the Neel 
order is developed which makes the concurrence becomes 
zero. The same treatment also is seen in the renormaliza- 
tion of entanglement entropy or concurrence. As stated 
in section (2) the renormalisation equations develops two 
different phases for the XXZ model which are separated 
by the unstable fixed point Ac = 1. For A > 1 region 
which falls into the Ising universality class, both the en- 
tanglement entropy and concurrence (Fig|2]) tends to zero 
as the iteration of RG equations increases, i.e the stable 
fixed point A — > oo is reached. 

A novel issue arises when the exponents of the scaling 
of Am and the minimum value of von-Numan entropy, 
as discussed in the previous paragraphs, are compared 
with the corresponding ones of the concurrence. The first 
derivative of concurrence tends to diverge at the quan- 
tum critical point. The position of minimum and the 
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FIG. 4: (color online) First derivative of entanglement en- 
tropy and its manifestation towards diverging as the number 
of RG iterations (steps) increases (Fig[3} • 



FIG. 6: (color online) The logarithm of the absolute value 
of minimum, ln(| dE/dA |a„), versus the logarithm of chain 
size, ln(A'^), which is linear and shows a scaling behavior. Each 
point corresponds to the minimum value of a single plot of 
FigH 
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FIG. 5: (color online) The scaling behavior of Am in terms 
of system size (A'') where Am is the position of minimum in 
FigH 



minimum value itself scale as = Ac + N~'^-^^ and 
I ^|a„ ~ iV° ''^, respectively, as the size of the sys- 
tem becomes large. It is clear that the exponents are 
almost the same for the two measures of entanglement. 
It indicates that the entanglement entropy and concur- 
rence scale in the same manner in the limit of large scale 
behavior and are associated to the long range behavior 
of the model close to the critical point. According to 
this intuitive picture of entanglement scaling and criti- 
cal properties of the model we would like to emphasize 



that the exponent is directly related to the correlation 
length exponent close to the critical point. The correla- 
tion length exponent, j/, gives the behavior of correlation 
length in the vicinity of Ac 

e - (A - Ac)-'^ (26) 

Under the RG transformations (Eq. p7|) ) the correlation 
length scales as ^ — > = S^/ub where ub stands for 
the number of sites in each block which is = 3 in our 
procedure. So, it is easy to look for the n*'* RG iteration 

- (A„ - Ac)-"- = iln% (27) 

which immediately leads to an expression for | 
in terms of v and ub 

l^k-iVV^ (28) 

The comparison with Eg. ((25)) demonstrates that the ex- 
ponent which governs the nonanalytic behavior of the 
entanglement entropy and concurrence in the vicinity of 
the critical point is nothing than the inverse of correla- 
tion length exponent. It should be noted that the scaling 
of the position of minimum, A„i (Figl5|) , also comes from 
the behavior of the correlation length near the critical 
point. As the critical point is approached and in the 
limit of large system sizes, not the thermodynamic limit, 
the correlation length almost covers the whole of system, 
i.e. ^ ^ N , and a simple comparison with Eq.f^H) results 
in the scaling of position as A^ = Ac -I- N^^^'^ . 
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V. FURTHER INSIGHT FROM QUANTUM 
GROUP 

As it has been carried out in Sec.jlll ini|) the (block) 
quantum renormahzation group is conceptually and tech- 
nically simple, but it may yield poor quantitative results. 
This was the main reason on the slow development of 
QRG through the 1980s in favor of powerfuU numerical 
Quantum Mont Carlo method. However there was a com- 
back to the RG in the 1990s as one of the most powerful 
method when dealing with the zero temperature proper- 
ties of the quantum systems. It was Wilson the first to 
associate the failure of the QRG to the role of the bound- 
ary conditions when applying to the tight bonding model. 
In fact the success of the density matrix renormahzation 
group (DMRG) as developed by White^ refers to the 
way in which it takes into account the effect of bound- 



ary conditions in terms of the quantum correlation in the 
ground sate of the system. Later, it was shown that the 
breakdown of the QRG is rooted on taking into account 
the entanglement in the ground state of the system^. 

However, Martin-Delgado, et aF" proposed a QRG 
prescription which implements the concept of the quan- 
tum group to the renormahzation group approach. In 
this approach the effect of boundary conditions have been 
imposed in terms of boundary magnetic fields on each 
block. The boundary fields cancel each other when col- 
lecting all blocks into the whole chain. This method was 
used to the describe the critical line of the XXZ model 
for < A < 1 and its basic idea relies on the "restoring" 
the rotational symmetry of the model by adding the ap- 
propriate boundary terms to the model. The open chain 
Hamiltonian is defined as 



JV 



1=1 



^q + q 



_z /q q \( „z _2 \ 



r 



Where q is an arbitrary complex number. 

To construct a real space RG, we have considered 
three-site block (Fig[T]) which is important in order to get 



a renormalized Hamiltonian of the same form as the orig- 
inal one. In this case the inter-block (H^^) and intra- 
block (H^) Hamiltonians are: 



J 



T -1 -1 

= I E + + 4/+l) 

1=1 



r 



In this way the correlations between blocks are taken 
into account in the coarse graining procedure in terms of 



the boundary fields of each block. The ground state of 
the block Hamiltonian is doubly degenerate and given by 



J 



IV'o) 



1 



--{^q'^'l TTi> + (g'/' + q-'^')\ TiT) - q-'/'\ iTT», 



V2(9 + g-i + l) 

1^^) = ^ ^ + (,1/2 + ,-1/2)1 ,-l/2| ^^))_ 



and the corresponding energy is cq = — 4(2 + g + q ^). The effective Hamiltonian with renormalized couplings is 
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7' 

H-ff = i. 

4 



N/3 , 



l) 



where 

, j'^ew . - +1 '''' 

The remarkable result of ea. (P5)) is that the coupling 
constant q when set as a pure phase, or alternatively 
A = — ) dose not flow under RG transformation 
which predicts correctly a line of critical models in the 
range |A| < 1, while J™, which is the value of J after m 
RG iterations, goes to zero in the limit where m ^ oo, 
which in turn defines the scale of energy. This means 
that in the critical region the entanglement measure, ei- 
ther concurrence or entanglement entropy do not evolve 
through the renormalization of the anisotropic coupling 
constant, i.e. it docs not rescalc the anisotropy param- 
eter (A' = A) since every point in this region is a fixed 
point. 

The density matrix is constructed from either lipo) or 
I'f/'o)- We then sum over site 1 and 3 degrees of freedom 
to get the reduced density matrix of site 2 and the rest of 
system. The von-Neumann entropy or the entanglement 
of site 2 and the whole of system is obtained by the eigen- 
values of the reduced density matrix (similar to Ea. ([M)) ). 
The entanglement of site 2 {Eg) has been plotted in Fig[7] 
versus the anisotropy parameter. For < A < 1 we ob- 
serve a decrease of the entanglement versus A which is 
the effect of anisotropy to reduce the quantum correla- 
tions whilt it is maximum at A = 0. Moreover, the 
reported values for < A < 1 do not evolve under RG 
transformation which shows to be the value as iV ^ oo. 
For A > 1 we get the evolution of Eg in terms of RG 
transformation. For higher RG iterations we get the zero 
value of entanglement which manifests the uncorrelated 
nature of a (classical) Ising state. The nonanalytic be- 
havior which manifest itself in the first derivative of the 
entanglement comes as the critical point Ac = 1 is ap- 
proached from the Ising (gapped) phase and is accompa- 
nied by a scaling with an exponents like what were seen 
before. 



VI. CONCLUSIONS 

In this work the relation between the entanglement 
and quantum phase transition via the renormalization 
group procedure is addressed. We have used the idea of 
renormalization group to study the quantum information 
properties of the XXZ chain. The entanglement entropy 
and concurrence as two measures of quantum correla- 
tions are used. In order to explore the critical behavior 
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FIG. 7; (color online) Quantum group manifestation of the 
evolution of Entanglement entropy at different RG iterations 
(steps) . 



of the XXZ model the evolution of both entanglement 
entropy and concurrence through the renormalization of 
the lattice were examined. As the number of RG itera- 
tions increases the entanglement entropy as well as con- 
currence develops two different values in both sides of 
the quantum critical point, where the phase transition 
between the spin fluid phase and the Ising like phase 
occurs. The phase transition becomes significant which 
shows a diverging behavior in the first derivative of the 
two measures. This divergence of entanglement measures 
accompanied by some scaling behavior near the critical 
point (as the size of the system becomes large). The scal- 
ing behavior characterizes how the critical point of the 
model is touched as the system size is increased. The 
fact that both measures of entanglement scale with the 
same exponent is a significant result of the large scale 
behavior of the model near the critical point. It is also 
shown that the nonanalytic behavior of both measures of 
entanglement is mirrored from the correlation length ex- 
ponent in the vicinity of the critical point. This inspires 
the behavior of the entanglement near the critical point 
is directly connected to the quantum critical properties 
of the model. 

Compensating the poor results of the QRG in the crit- 
ical region of the XXZ model, we resort to the quantum 
group method in order to take into account the corre- 
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lation between blocks. In this way the massless region 
is truly captured and a realistic feature of the evolution 
of the concurrence or entanglement entropy is inspired. 
However, the nonanalytic behavior and the scaling be- 
havior is the same as before as the critical point of the 
model is approached from the gapped phase. 

The approach presented here can be used to get the 
quantum information properties of quantum systems in 
an easy way. In other words we get the properties of 
a large system dealing with a small block which make 
it possible to get analytic results. However, the quan- 
tum renormalization group usually suffers from the poor 
quantitative results. One should get more accurate 
results using the proposed idea with a more complex 



numerical method like density matrix renormalization 
group or exact diagonalzation method. 
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